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Ionograms  recorded  at  Brighton,  Colorado  during  the  last  sunspot  maximum  are  still  being  studied.  Samples  show 
that  the  high  quality  and  high  temporal  resolution  of  the  data  may  yield  new  information  about  the  dynamics  of  the 
ionosphere,  especially  in  the  F-region.  F-region  parameters  from  these  data  show  that  this  part  of  the  ionosphere  is  highly 
variable,  with  short  periods  of  the  order  of  fractions  of  an  hour. 


This  report  explains  why  the  maximum  usable  frequency  (MUF)  (3000)  appears  to  be  the  most  appropriate  indicator 
for  F-region  variations. 
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INTRODUCTION 


During  the  last  sunspot  maximum  in  1980  and  1981,  ionograms  were  recorded  at 
Brighton.  Colorado  with  a  new  digital  ionosonde  developed  by  the  Space  Environmental 
Laboratory  of  National  Oceanic  and  Atmospheric  Administration  (NOAA).  cosponsored  by 
the  National  Science  Foundation.  At  the  time  the  program  was  terminated  in  1981.  onK  a 
small  percentage  of  the  data  recorded  had  been  analyzed.  The  results  of  those  samples 
indicated  that  the  high  quality  and  high  temporal  resolution  of  the  data  may  yield  new 
information  about  the  dynamic  of  the  ionosphere,  especially  in  the  F-region.  This  justified  a 
more  systematic  and  detailed  studs  of  all  the  data  available.  All  the  F-region  parameters 
derived  from  the  data  show  clearly  that  this  part  of  the  ionosphere  is  highly  variable,  with  short 
periods  ot  the  order  of  tractions  ot  an  hour.  T  hese  observations  raise  questions  about  the 
physical  nature  ot  the  phenomena  like  cause,  origin,  propagation  speed,  and  direction.  I  hev 
also  provide  some  indirect  information  about  the  effects  on  high  frequency  (HI)  propagation 
due  to  variable  tilts  ot  contours  ot  constant  electron  densities.  Mainly  tor  application  in  the 
H 1 -propagation  area,  but  also  tor  correlations  with  other  geophysical  parameters,  it  is 
desirable  to  have  a  parameter  that  can  give  a  first  order  estimate  oi  the  F-region  variability.  As 
will  be  discussed  later,  the  maximum  usable  frequency  ( Ml  Fit 3000)  appears  to  he  the  most 
appropriate  indicator  for  the  F-region  variability. 


We  will  show  that  F-region  oscillations  are  present  all  the  time  with  variable  intensity, 
that  those  variations  are  well  represented  by  the  Ml  Ft 3000).  and  that  correct  estimates  ol  the 
rate  ot  change  of  the  Ml  F( 3000)  and  therefore  of  the  F-region  variability  require  a  sampling 
rate  of  at  least  10  to  12  ionograms  per  hour,  much  faster  than  the  standard  rate  of  4  ionograms 
per  hour.  Since  very  few  ionogram  sequences  with  a  high  sampling  rate  exist,  we  show  plots  of 
the  MI  F(3000)  for  all  the  rapid  sequences  taken  during  the  lifetime  of  the  program  in  the 
appendix  A.  For  comparison  we  also  show  most  of  the  1 5  minute  interval  data  m  appendix  B 


THE  DAT  A 


A  short  outline  of  some  of  the  features  of  the  ionosonde  is  appropriate  lor  a  better 
understanding  of  the  quality  and  the  information  content  ol  the  data.  A  more  complete 
technical  description  was  given  bv  (irubb  ( I9’9(  lour  receiving  antennas  located  in  the  corners 
ot  a  square  with  a  diagonal  distance  ot  100  m  were  an  essential  feature  of  the  system.  Felloes 
reflected  from  the  ionosphere  were  recorded  at  two  antennas  simultaneously  by  a  two-channel 
receiver  capable  of  recording  amplitude  and  phase  of  the  signal.  All  data  wc  ■  digitized  ami 
then  recorded  on  magnetic  tape.  A  basic  data  set  was  obtained  by  sounding  a  pair  of  fre¬ 
quencies  (differing  by  8  kHz)  and  recording  the  echoes  at  all  four  antennas,  which  required  the 
transmission  of  four  consecutive  pulses.  T  he  switching  of  antenna  pairs  and  frequencies  was 
performed  in  such  a  way  that  the  variation  of  the  phase  with  frequence  could  be  separated 
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Usually  the  time  inters al  between  consecutive  pulses  was  20  msec.  This  configuration  in 
space,  frequency  and  time,  and  the  high  accuracy  ot  the  phase  of  1  degree  or  better  permits  the 
evaluation  of  the  following  quantities  (Paul  et  al..  1974): 

the  angle  of  arrival  of  the  echo  from  the  variation  of  the  phase  with  antenna 
location, 

the  Doppler  frequency  or  Doppler  velocity  from  the  change  of  the  phase  with 
time. 

improved  estimate  of  the  virtual  height  from  the  change  of  phase  with  frequency. 

The  ionosonde  operates  on  a  quasi-logarithmic  frequency  scale  and  depending  on  the 
frequency  density  selected,  a  high  quality  ionogram  can  be  recorded  in  less  than  a  minute. 

During  the  period  of  time  when  the  data  in  question  were  recorded,  the  digital 
ionosonde  was  in  its  final  stage  of  development  and  extensive  testing  of  hardware  and  software 
took  place.  For  this  reason,  data  were  recorded  in  an  irregular  fashion.  For  example,  very  few 
data  were  recorded  during  the  summer  months  of  1980.  Also,  during  this  time  it  became  more 
and  more  evident  that  the  standard  sequence  of  recording  ionograms  in  15-minute  intervals 
actually  undersamples  the  more  or  less  regular  short-term  variations  of  the  ionosphere. 

1  hercfore.  the  majority  of  the  rapid  sequences  of  ionograms  were  taken  toward  the  end  of  the 
progra  m. 

I  he  standard  automatic  processing  routine  used  for  this  study  provides  the  following 
parameters  for  each  ionogram: 

foF2.  the  critical  (penetration)  frequency  of  the  F-layer  for  the  ordinary  component 

fxF2.  the  critical  (penetration)  frequency  of  the  F-layer  for  the  extraordinary 

component 

hmoF2.  estimate  of  the  height  of  the  F-region  maximum  derived  from  the  ordinary  echo 
trace 

hm.xF2.  estimate  of  the  height  of  the  F-region  maximum  derived  from  the  extraordinary 
echo  trace 

YmoF2.  estimate  of  the  half  thickness  of  the  F2-layer  derived  from  the  ordinary  echo  trace 

Ym.xF-2.  estimate  of  the  half  thickness  of  the  F2-layer  derived  from  the  extraordinary  echo 

trace 

Ml  f-( 3000).  (the  maximum  usable  frequency  over  a  7()()()-km  path),  derived  from  the  ordinary 
echo  trace 

IMl'F.  the  frequency  of  the  tangential  point  of  echo  trace  and  MUF  curve 
hMl'F.  the  virtual  height  of  the  tangential  point  of  echo  trace  and  Ml'F  curve 

/I  the  zenith  angle  averaged  over  the  ordinarv  echoes  at  six  frequencies  in  the 

vicinity  ot  the  tangential  point 

\X.  the  azimuth  angle  averaged  over  the  ordinarv  echoes  at  six  frequencies  in  the 

vicinity  of  the  tangential  point 

the  highest  frequenev  where  an  F-region  echo  is  observed 


HI  s. 


While  all  the  F-rcgion  parameters  listed  above  do  reflect  the  changes  in  the  F-region. 
most  of  them  do  not  qualify  as  a  good  variability  parameter  for  various  reasons.  For  example, 
the  height  of  maximum  and  even  more,  the  half-thickness  are  both  extremely  sensitive  to  small 
deformations  of  the  layer  and  temporal  changes  of  those  parameters  cannot  be  easily 
interpreted  as  changes  in  height  or  electron  density.  This  is  also  true  for  the  quantities  fMl'F 
and  hMUF.  For  these  reasons  the  following  discussion  of  a  variability  parameter  will  be 
limited  to  the  parameters  foF2  and  M  L  F( 3000). 

Approximately  23.000  ionograms  recorded  at  Brighton  were  processed  to  obtain  for 
each  one  all  the  parameters  listed  above.  The  same  type  of  digital  ionosonde  became 
operational  in  late  1980  at  White  Sands.  NM.  Unfortunately,  most  of  the  data  recorded  there 
are.  for  various  reasons,  not  well  suited  for  automatic  computer  processing.  This  is  especially 
true  for  the  overlapping  period  of  time  w  hen  both  systems  were  operational  and  only  very  few 
data  were  available  for  comparison  of  the  two  sites.  For  these  reasons  no  White  Sands  results 
will  be  reported  here. 

THE  Ml  F(3000) 

The  maximum  usable  frequency  over  a  3000-km  path  can  be  derived  from  a  vertical 
ionogram.  The  process  is  based  on  the  relationship  between  the  vertical  and  oblique  frequency 
to  be  reflected  at  a  given  height  for  a  given  electron  density  distribution.  If  the  electron  density 
is  a  function  of  height  only,  it  is  sufficient  to  know  the  vertical  ionogram  to  establish  this 
relation  (Smith.  1939).  In  the  traditional  analog  process,  the  MUF(3000)  is  determined  by 
overlaying  a  set  of  transmission  curves  on  the  ionogram.  The  virtual  height-frequency  relation 
of  those  curves  was  published  by  URSI  (Union  Radio  Scientifique  Internationale)  in  the  form 
of  a  table  (Piggott  and  Rawer.  1972).  An  empirical  algebraic  expression  for  those  curves  of  the 
form 

r  =  MUF  f  =  (a  -  b'h'ls/h7 

was  derived  by  Paul  ( 1984).  Here  h'  is  the  virtual  height. 

a  =  67.629  and  b  =  -0.0148 

For  a  given  ionogram  h'(f)  the  MUF(3000)  can  then  be  determined  in  a  least  square  fit 
using  the  above  expression.  The  coordinates  fMUP  and  hMUF  of  the  tangential  point  between 
the  echo  trace  and  the  MUF  curve  are  also  obtained  in  this  process.  We  found  that  for  most 
ionograms  the  frequency.  fMUF  is  on  the  average  approximately  1 1  r<  less  than  the  critical 
frequency.  foF2.  Assuming  a  parabolic  model  for  the  F-layer,  this  frequency  ratio  corresponds 
to  a  true  height  somewhat  above  the  middle  of  the  lower  half  of  the  layer,  or  about  45r(  of  a 
half-thickness  below  the  height  of  maximum.  This  relationship  gives  a  strong  argument  for 
considering  the  variations  of  the  Ml  F( 30()0)  to  be  representative  of  temporal  changes  of  the 
lower  half  of  the  F-layer.  An  alternative  parameter  representative  of  variations  of  the  F-layer  is 
its  maximum  electron  density  or  critical  frequency.  foF2.  In  our  experience,  three  factors  favor 
the  M  U F( 3000 )  over  foF2.  The  critical  frequency.  foF2.  is  obtained  by  extrapolation  based  on 
some  model  assumption  and  is  very  sensitive  to  small  changes  in  the  shape  of  the  actual  layer. 

I  his  means  that  apparent  variations  in  foF2  are  likely  to  be  more  influenced  by  changes  of  the 
'lope  of  the  profile  in  the  range  used  for  the  extrapolation  rather  than  by  changes  of  the 
maximum  electron  densitv  Anothei  factor  in  favor  of  the  MUF(3000)  is  the  observation  that 
in  the  presence  of  weak  spread-F  the  echo  trace  in  the  vicinity  of  foF2  is  often  poorly  defined. 


I  NT  HRPRHTATION  OF  THF  Ml  F( 3000)  VARIATIONS 


Variations  ot  the  Ml  1(3000)  can  he  caused  by  changes  in  height  or  electron  density  ot 
the  F-layer  or  a  combination  ot  both.  In  general,  we  will  not  attempt  to  distinguish  between 
those  effects,  but  we  will  study  a  few  examples  in  greater  detail.  Figure  3  shows  an  example  of 
M  IF  variations  derived  from  ionograms  recorded  with  different  sampling  rates,  first  for 
almost  3  hours  at  a  rate  of  30  ionograms  per  hour;  later,  only  4  ionograms  per  hour  were 
recorded.  \  comparison  of  the  two  data  sets  gives  a  strong  indication  that  the  lower  sampling 
rate  definitely  undersamples  the  true  variation  of  the  F-region.  The  oscillations  in  the  morning 
have  periods  of  approximately  20  minutes  with  amplitudes  in  the  order  of  1  MH/.  While  the 
length  of  the  periods  show  n  is  quite  common  for  all  the  data  recorded  and  analy  zed,  the 
amplitudes  are  often  much  larger.  This  is  shown  in  a  different  way  in  figure  4  where  large 
changes  over  b-nnnute  intervals  of  the  entire  F2-laver  echo  trace  become  clearly  visible.  1  he 
Ml  1  values  for  the  sequence  were  22.33.  21.35.  and  22.09  MH/  respectively;  first,  a  decrease 
(it  I  MH/  and  then,  an  increase  of  OT  MH/  over  the  short  interval  of  6  minutes  each. 
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A  closer  look  at  the  upper  curves  in  this  figure  reveal  a  time  lag  between  the  variations 
at  the  higher  and  lower  frequencies.  For  example,  the  last  minimum  between  1500  and  1000 
appears  later  at  the  lower  frequencies  than  at  the  higher  ones.  The  same  is  true  for  the  pre¬ 
ceding  maximum.  We  also  see  that  relatively  short  periods  of  approximately  12  minutes  are 
present  between  1430  and  1500  with  larger  amplitudes  at  higher  frequencies  than  at  lower 
frequencies.  Those  short  periods  are  also  clearlv  visible  in  the  MTF;  curve. 

The  high  accuracy  of  the  MIT  (3000)  permits  the  estimate  of  its  temporal  derivatives 
by  simply  taking  the  divided  differences  of  consecutive  data.  1  wo  examples  of  this  quantity  are 
shown  in  figures  6  and  7.  The  top  portion  of  figure  6  shows  the  MTF(3000)  from  noon  to 
midnight  for  12  March  1981.  the  temporal  derivative  of  the  Ml'F(300())  is  plotted  in  the  middle 
of  the  figure  The  ionograms  were  taken  in  3-minute  intervals.  The  MTT  curve  shows  a 
relatively  smooth  behavior  with  some  oscillations  of  various  durations  superimposed.  I  he 
variations  become  much  more  visible  in  the  temporal  derivative.  We  notice  that  many  times  the 
magnitude  of  the  derivative  can  be  comparable  or  even  larger  than  the  average  rate  of  decrease 
during  sunset.  In  order  to  demonstrate  the  aliasing  effect,  we  eliminated  from  the  original  set 
all  the  data  except  those  recorded  at  the  full  15  minutes.  The  temporal  derivative  derived  from 
those  data  is  shown  in  the  bottom  portion  of  figure  6.  The  effect  of  undersampling  is  visible  in 
different  ways.  For  example,  the  half-hour  oscillations  between  1200  and  1400  are  still  present 
in  the  15-minute  data,  but  their  amplitudes  have  only  approximately  halt  the  value  of  the  origi¬ 
nal  data.  Some  rapid  changes  of  the  derivative  take  place  between  2000  and  2100  as  seen  in  the 
original  data.  This  information,  however,  is  lost  completely  in  the  15-minute  data.  Figure  7 
shows  the  continuation  of  this  data  set  and  includes  sunrise.  Again,  very  large  magnitudes  of 
temporal  changes  are  clearly  visible,  some  even  larger  after  sunrise  than  during  sunrise.  A  com¬ 
parison  of  the  time  derivatives  of  the  original  data  with  the  15-minute  data  e  g.,  between 
0700  and  0900  -  confirms  the  earlier  finding  that  significant  information  can  be  lost  by 
undersampling. 

Angle  of  arrival  measurements  provide  some  evidence  indicating  that  the  temporal 
variations  observed  are  not  simply  upward  and  dow  nward  movements  of  the  ionosphere  as  a 
whole,  but  rather  propagating  weak  deformations  of  the  F-region.  most  likely  acoustic  gravity 
waves.  If  we  assume  a  propagation  velocity  comparable  to  the  speed  of  sound,  we  can  translate 
the  temporal  variations  into  an  estimate  of  the  horizontal  gradient  of  the  MTF(3000).  On  the 
other  hand,  it  turns  out  that  for  a  given  location  at  midlatitudes  the  relative  velocity  of  the 
day  light  boundary  is  also  close  to  the  sound  velocity  depending  on  the  season  of  the  year. 
Therefore,  by  comparing  the  temporal  variations  at  a  given  time  of  the  day  with  the  temporal 
variations  during  sunrise  or  sunset  we  can  icasonably  assume  that  this  ratio  is  also  in  first  order 
valid  for  the  corresponding  horizontal  gradients,  where  the  ones  for  sunrise  and  sunset  are 
relatively  well  known. 

Tables  la  and  lb  give  some  statistics  for  all  the  rapid  sequence  ionograms  ( 10  or  more 
ionograms  per  hour)  that  were  recorded  during  the  duration  of  the  program.  We  computed 
lower  quartile  ( l.Q).  median  (MFD),  and  upper  quartile  (I  Q)  of  the  absolute  value  of  the  time 
derivative  of  the  MI  F(3000)  for  three  time  intervals  of  the  day  :  sunrise,  daytime,  and  sunset 
Nighttime  values  are  not  shown,  since  frequently  spread-F  condition  caused  relatively  large 
uncertainties  of  the  MUF(3000)  values,  and  separation  of  uncertainties  from  actual  variations 
would  have  been  difficult  and  time  consuming.  For  nights,  when  no  spread-F  was  present,  the 
data  show  that  the  variability  at  night  can  be  as  high  as  during  the  day  For  a  given  day.  the 
median  values  are  usually  largest  for  sunrise.  1  here  is  a  definite  asymmetry  between  sunrise  and 
sunset.  We  also  see  that  median  values  tor  daytime  are  not  much  smaller  than  those  tor  sunrise 
and  that  the  day  time  upper  quartile  is.  on  the  average,  close  to  the  sunrise  median.  Similarly, 
the  sunset  upper  quartiles  are  not  much  different  from  the  daytime  upper  quartiles. 


I  ublc  2.  based  on  15-minute  data,  has  been  added  for  comparison  only  and  the  data 
should  be  used  with  caution  As  we  have  shown  above,  we  have  reason  to  believe  that  the 
relatively  small  day  time  medians  are  the  result  of  undersampling  and  do  not  necessarily  reflect 
lower  gr  avity  wave  activity.  In  addition,  it  should  be  mentioned  that  typically  the  sunrise 
period  lasts  for  about  2  hours  and  with  15-minute  sampling,  only  9  data  are  available,  which 
may  not  be  a  large  enough  sample  for  a  good  estimate  of  the  statistics  listed.  Generally,  the 
data  in  this  table  give  the  long-term  trend  (hours)  of  the  Ml  T  ( 3000)  and  are  less  influenced  by 
t he  shorter  v anationv 


I  able  la  Absolute  values  of  temporal  changes  of  the  MIT(3()00)  in  MH/  hour  for  all 
rapid  sequence  lonograms  from  i.5  F  ebruary  1980  through  3  January  1981. 
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I  a h i c  I  b  Vhsolute  value'  ot  temporal  changes  ol  the  M  l  H l(H)O)  in  M  H/  hour  la  all 
rapid  sequence  lonograms  trom  .lainiarv  I'JM  through  16  September  I9KI 


(  ONC  I  I  SIONS 


We  consider  three  results  ol  this  studs  as  highly  significant: 

(a!  Oscillations  ot  the  F-laycr  are  taking  place  more  or  less  all  the  time  with  varying 
intensity  and  they  are  ssell  represented  by  the  temporal  variation  ot  the 
Ml  1(3000). 

(hi  Since  the  periods  ol  those  oscillations  are  relative!)  shott.  correct  estimates  ot  the 
temporal  variations  ot  the  1  -region  parameters  can  only  he  obtained,  il  at  least 
10  to  12  lonograms  per  hour  are  recorded. 

(c)  lemporal  variations  and  associated  tilts,  comparable  with  those  during  sunrise 
and  sunset,  can  occur  at  any  time  ot  the  das  or  night. 

1  or  further  discussions  we  assume  a  propagation  velocity  ot  300  m  s  and  a  period  ol 
20  minutes  I  his  would  then  correspond  to  a  wavelength  ot  360  km.  il  we  can  assume  a 
periodic  spatial  structure.  Consequently,  spatial  sampling  should  take  place  with  grid  distances 
ol  less  than  halt  a  wavelength  or  less  than  1X0  km  il  a  true  picture  ol  the  ionospheric  structure 
is  to  be  obtained.  With  the  same  velocity  in  a  linear  approximation  a  rate  ot  change  ot  the 
Ml  1(3000)  o|  10  MH/  hour  would  correspond  to  a  change  in  a  horizontal  direction  ol 
approximately  I  MHz  per  100  km.  As  mentioned  earlier  the  quasi-periodic  changes  ol  the 
Ml  1(3000)  are  mainly  caused  by  changes  in  the  virtual  heights  and  to  a  lesser  degree,  by 
changes  in  the  electron  density.  I  he  variations  ol  the  virtual  heights  correspond  to  changes  ot 
the  true  reflection  height  over  an  area  with  dimensions  ot  the  order  ol  a  traction  ot  a  wave¬ 
length.  By  this  relation,  the  temporal  variations  ol  the  Ml >(3000)  are  also  indicators  lor  tilts. 
In  figure  2.  the  fourth  curve  from  the  bottom  follows  the  virtual  height  very  closely  where  the 
\1l>  is  determined  (hMl'f  ).  Comparing  this  curve  with  the  Mil  curve  in  the  middle  section 
ol  the  figure,  we  find  that  an  amplitude  of  the  Ml  F( 3000)  oscillation  of  approximately  0.5 
MHz  corresponds  to  a  variation  of  the  virtual  height  with  an  amplitude  of  approximately  15 
km  Hy  a  conservative  estimate  a  virtual  height  interval  of  15  km  is  equivalent  to  a  true  height 
interval  of  approximately  5  to  ?  km.  Such  a  change  ol  the  true  height  over  a  quarter  of  a 
wavelength  would  then  represent  an  average  tilt  lor  a  surface  of  constant  electron  density  ol 
3  4  degrees.  Angle  off  arrival  measurements  verity  that  the  equivalent  angles  ol  deviation 

from  vertical  propagation  are  observed  frequently. 

While  this  type  ol  consideration  can  give  some  lirst  order  estimates  ol  some  quantities 
defining  the  temporal  and  spatial  structures,  other  very  important  parameters,  such  as  the 
propagation  velocity  of  an  acoustic  gravity  wave  in  its  magnitude  and  direction,  cannot  be 
obtained  by  single  site  vertical  sounding  knowledge  ot  this  velocity  could  have  important 
applications  ranging  from  locating  the  sources  of  gravity  waves  to  directional  updating  ol  Hf  - 
propagation  conditions.  I  he  Doppler  velocity,  which  can  be  observed  by  an  lonosonde. 
represents  only  the  radial  component  ot  the  temporal  change  of  the  phase  path  It  lags  in  time 
behind  the  group  path  velocity  (Bennet  and  Dyson.  19X6)  and  for  these  reasons,  cannot  be 
easily  interpreted  as  a  propagation  velocitv  of  gravity  waves 

It  has  to  he  pointed  out  that  the  data  set  available  lor  this  study  was  insufficient  to 
come  to  definite  conclusions  about  the  seasonal  or  even  the  solar  activity  oi  global  dependence 
of  the  F-region  variability.  The  data  show,  however,  that  the  F-region  was  highly  variable 
during  the  recording  period  and  several  parameters  showed  olten  shorter  periods  and  larger 
amplitudes  than  previously  anticipated  We  further  come  to  the  conclusion,  that  no  database 
exists  to  describe  the  short-term  dynamics  of  the  l  -region  and  its  corresponding  medium  scale 
structure  on  a  larger  scale  in  space  and  time  Answers  to  many  remaining  questions  have  to 
come  from  a  continuous  observational  program  with  sufficient  temporal  (.it  least  10  lonograms 
per  hour)  and  spatial  (distances  not  more  than  =  200  km)  resolution 


C  OM  MEN  I  S  ON  APPFN DIC'FS  A  AND  B 


Appendices  \  and  B  contain  plots  of  the  Ml  F(3()00)  tor  most  of  the  data  available. 
Appendix  A  shows  the  data  tor  the  rapid  sequences  ot  ionograms  (8  or  more  ionograms  per 
hour).  Appendix  B  shows  the  standard  sequences  (less  than  8  ionograms  per  hour).  Some  of 
the  data  shown  appear  unreasonably  large  or  small.  In  those  eases,  an  accurate  determination 
ot  the  \alue  was  not  possible  due  to  spread-F  condition  or  multiple  reflections  due  to  intense 
deformation  ot  the  laser.  The  continuity  of  the  data  derived  from  the  rapid  sequences  of 
ionograms.  in  most  cases,  gives  a  dear  indication  if  deviations  from  the  mean  trend  represent 
real  variations  or  poor  data.  This  distinction  is  often  not  possible  for  the  standard  sequence 
data.  1  he  Ml  F( 3000)  values  are  set  to  zero,  if  for  technical  (e  g.,  incomplete  ionogram)  or 
natural  reasons  (blanketing  sporadic  F)  a  correct  value  could  not  be  obtained. 

It  has  to  be  pointed  out  that  each  Ml  F  value  shown  was  derived  from  an  individual 
ionogram  and  no  smoothing  or  filtering  was  applied. 
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